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1 . 

£^**1 ) has . t t*d a n*rlcal f- ' *t \*\p* *U £t*?***\*l 

equation* 4*1##* U*es rational Pad* C^rVeTf** ts to ipr^HU fvwUfM 
of the f&rwmrd shift operator S to generate rao«*si«r> taraa&aa,, fcw hmm 
deess*strated tha merits of using such foml&s In solrin# first e rt^ 
•quatiiaij and has exhibited certain stability problems related ^»nu ( 
Tha purposa of this pap ar is to extend tha study of this me%hma m 
tha spaoial oategoiy of seoond ordar differential aquations in vtUb she 
first derivative doas not appaa r. A computer program is darali^i 
which ganaratas a nt=£er of raoursion formulas, applias than nMiaiTa. 
ly to solve tha initial value problan y«+ y * 0, y((} * 1, y»(0) * 0 f 
and producas graphs of solution arror vs number of iterations, Tha ba~ 
tic program is than modified to solve tha Coulonb Wara Equation *j ' + 

^ ):o with initial conditions appropriata to tha sal ao tad 
values of tha parameters y and L, Tha Runge-Kutta method is uaad to 
compute eight values of tha independent variable to start tha solution, 
finally, tha results of these computations are discussed in detail, 

Sone of tha advantages of this method over others in common use ares 

a, A rational convergent is generally a batter approximation 
to a function than is the corresponding truncated series. This will be 
discussed further in the following section, 

b. Only about one half the total manber of tabular differ- 
ences required by methods using polynomial approximation a are needed 
in this emthod for comparable accuracy. This is because differences 
ar# distributed over both sides of the recursion formyla, 

o. Since differencing is divided between both sides of the 
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*«U±-r «|4» Ift 'Itt*. 



•tu* 



nw~ T c/ Mqdryi «*4* »-"■ * |t — tlj, O* (?•# «*T1 

4. Wtf€i« Costal** IMTIMT ?■*• ■ 90 m 

u--~«Ai »* »* W? »«rvtt« > • dJ cl 1*1 - r — E>« i *>*J 

AAJU r » * *■ U Vorign *. 

2. I«TlT«CW, FBP wtlM . 

*<f'p*«* it 1* d*eJp»d to — W* ’»*• kAltiftl t«Iv« ^n»Ln 

/*(*) • f^r»(x)), y(x.) - jr. , jr'U.) « y; . 

I#r fonuuUy .it mca «Ai# 

\ yO it *t 

w . nrt c tb« l*ft «ld* 

[*../•• io|-l ' •)!* •*•>- ^ U ' * 

VI 

>.i_— ^ T r-Vn C * ?bcr<9 of variable 

t ♦ * + te ]t d f * hdt x 

t • c rf t .Si , dt * fcrft, 



*id da! ini/*< ->• focv*rd ihift operator E irJ^n that I f(x) * f(i 
» pi t_Afc rifQt (traatinf 1 a# an alftibralo ^feol)7 








A » ,* - 



J J 
■ • ^ 




• •• lA* r- 1*1 *** tm 

V ,-'••< 

S.U1M1.1 1- A '\ '' 

u-oj. a,: .ilf. , . i. . .Cii. , 

J |il- ,t/ "if 1 \ »'5*J #-• 7 *.,# 

In t«i*a ot the backward dlffarano* oparat or », d»^U« ( t^r ▼ (*,■ • 

y, - y,-., » is 



H k If 4-t-ti^s' -_i_ 

. „ Ing through by 

'V 5 {*vfcr] t. * l ' (T -' , ’V * i^rs ■ 



Multiplying through by e 
V 

since ? operating on a constant yields aero. 

For simplicity of notation we let i x * v*yV «w* write 




**> 

X V ^ ' f 

Let a # + a t T + a^T + ... be the power series expan* lee* of* . 



The power series can be approximated by a rational free tip* b> 

the following technique, a + a t T + a.v x + ••• + a_V£ ik * 8 * ■ r * VI 

Multiplying each side by the denominator of the above f r**~tirr» iq) 

equating coefficients of the like power terms in v we fmmmkm jH-q+1 

linear equations In the p+q+1 unknowns b # ... ^ , o • . .o f . 

b # - a, 

b, a a, + a, o, 

b,« a*+ a, c + a.Oj. 

b i* a,+ ♦ a, c 4 + a # c i 

The solution of this © stem gives us the o*a£flalwts of 
*1 fraction which agrees with the series o*-t to P+qpM 
It should be r*oted ffeat a given po f *rmsalal of 'err** r h* b< ■**»«» Lp-* 
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to of r + i dj*Uwt fr «rw ♦ • , 1 . - r B ,,_ 

w-dul *fUr 1, 2, ...» n t*™*» :S*“ » 

of polynomial fraction* wtu** »*n »♦ *rr«mgi 1 a* * t *d* Uel* •* *•«-<-.- 

, >» .i ' w< i 4 

*1 conv*rg«vti. Fig. 1 la th, Pad* t*bl* asa?Qlai«j wi*i» f — 
wbara a third dagraa polynomial la uaad to f«^«rat« t>i latliMl 



oonrarger.ta. 

Tha power series expansion 



-i^V . 



tens, Is the sixth degree polynt^dal 

1 \l ^ W |V& 11 5** 

The coefficients of this polynomial were used as tm 

/ 

S UCSD PADS (Naval Postgraduate School Library Routine) to dbhal* 



(5) 



coefficients of the 26 Pade oonver gents whose numerates rnmt W% 

are polynomials of degree q and p, respectively, and Mere p*q 2 6. 

Pig. 2 lists these coefficients, rounded to three signiflMIt fif9ee. 

As stated above, each of these convergent s ie identically m\ml to 
(3) out to p+q+1 terms. Moreover, a given convergent expanded alpv^veje^ 
ally into a power series generally continues to agree quite well wi*k 
the related infinite series for a n«ber of terms beyoed the (ptq+1)* 1 ! 
i.e., a rational convergent generally is a better approximation to a 
function than is the truncated series from which the omr.ver^t ms 
derived. 

j ityi -») ] 1 

Desiccate a partioular otmwerfwat to ( ^ by 

ftNf-L. AUiAH. * ftAL , .V.2 1 ■■» V r,Jii *1. . (*) 

a ♦ a »n> v a 4 & j|tr 1 f * 

Suppose re wish to find an apprettMte numerical nl<aUx) ef the aqwsfc3 ■» 

r j.«i ihr) } x 

/" » f(x,/), fchstitatlfcA (4) for ^ — j In eq-ction (2) we Mte 

i * * H»UxlT. * ) . 

/ \ *t>< . i ♦ A* i* 

n IMi * , a , » • • • * V - ' 1 - 



«• 






l>* 






<n w. 




* >>jq* - •r * 

... w- (i ■ ' <- '' ■ I 4 -. -‘ fv- 

.«• . ,*'[<£• v. -|' 4. '*> **'••"'• • • (ii 

-Wk-'U **e ■ ->Tj v. O) »* 

QU> »JE. 

ro) i\*oU’D 

44jt*ti t iLini ts* in (5) *• 

<j(u)r, - | ♦ . ,.: * h ^P(o)t, -P(i)C,. ♦ ...: . *** 

t % - ••* • ®(j|)%.- f + h*£p( 0 )f, * •••- p (p) f ' r] ^ • 

{]ifi | ^ I r In » «~» * of OOfiffidWlts of ^ ’• t nr*»«.** l Mi > . 

of previously oMputed values of i and f. R*calli«« that i » v y, - 
y> _ 2y, + y, , we My now calculate y ( • *,+ 3y,., - 7, x . The value of 
in (6) will have b*wn ooeseited fro* a flret appr^S^sti^ to 
y f«*»*rat«r» by a predictor forwula. 



(*) 



3. stability C neidaratlons . 

The recursion foreula jenareted froc a given Pad* ewer* ««It is *S 
tne forst Q(0) t, • Q(l)r<.. - . . . i Q(q)*‘^ + h » ... /P(p^T, (7) 

whet* Q(i) and P(i) are as defined in Section 2. 
tklbetittfting v l y, » y, - 2| .,+ 7i j. fori, in (7) we obtain 

q(o) (y, -5? , , ) » QU) (y«* -2y, ,+y, ,) - ...^Q(q) 

-2y. f +y^.,) * h^p(o)f t - p(i)f ,+ ...i p(p)t,-t‘] 

Udlr4^iyin< .,<^1 collectius te«-*e, this U«su > 

QL 43 )f, • [q( 1 ) + Si( 0 )^ y, - jo( 2 ) + *U) + Q(°)J J,. a + *** 

* (-1 )**’ | Q(q) + A(%-1) ♦ 3(q>2)] 7^ ♦ (--)* [**(*) ♦ a ^ l ^j I'-f t 

+ (-1)* 4Vh)» t ^ i |**( c )A - • ••♦ (- 1 ) p (p) f ~r ] 

Tbs *s were} elUtiun of la*. * it ?f*w» wi«*U»* is ot tbs f« 
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j . - : f f , , r , • ■ > f - 

rft. 1 * of the la — Mlm.it*. Win 

later) .<si the c, ere 

CmaKiar the 2ad order Pade odtfH»t«- ,1 „ ^* B "** # 15 

Ui fl|, *). Frn» the discussion of 9rati«*» 2 *» teoet U>*1 t*e *wsa«t^ 
f omul a gmerated by this coTTergwt will have the foU»-<to{ 

Q(0) « It Q(i) - 0, i * 1 

p(o) *4? ■ -!{* P(2) P(i) “ °t 1 * 3 

general recursion formula (8) then becomes 

y, ‘ 2 y,„- y,.*+ • 

In the Initial valua problem (1) it la seen that f , ■ -y, . 
Substituting in (9) end o oil eo ting terms} 

(i +4 hl )y; * U-gJ? )7,j - (l+^h^y, * 

or. 

(12+h 4 )y-- (24-1 Oh* )y,_ + (12+h" )y,^« 0 

t ml »)* 



(J > 



finding roots: 



/v = 



k( iHh*) 



fll-5^ t a/ V 144- 3-4 1,* 
/2LV-A*- 



Ignoring tans of the order of h we have 

/v 2; / 2.-^ <3 f ^ i / 

»a- 



( 10 ) 



(ID 



The general solution of (9) is of the fora 
o. r> o,r' 

deteteme the oieffioisnts o,‘ froa '-he eolation y • we x *f ' 
Initial ralue problem. 



S, * cr>o.r ( < 1 



(initial c3-ii*-i;y*e f.’ 



y, • o r. + e, r 



a a ^ 

1 - — (y • «>t(i • 1 -?+*[••• . tm eaU 

3 ralaee *T a «oad 



ft 










4 1 



J* 



&«,« we !*«v«S 




(*£• 



ana can b* writt«fi In trtfwowtrlo ror* t 





wh«r* 



& z OJ^*~ 





o 



( 13 ) 



Cot «*<4Jt h, h^in afalaH rtg« 

^ - (V <♦*.») ^co»-*>4 - -^- S 

and this solute on cocperfes to the exact solution. 

It can bo from the grafts that Formulas 10, 11, 12, 15 , 2?, 
ec«l 25 aro tho moot pM sing fry* tho standpoint of orror 

«-><*Ltfication. Onfortaaately, their associated characteristic eq-*tio** 
\Lro< tevenife dicrse) are too intractable to pe^it analysis 
It wmild e quit* alible, loirever, to rbtain the roots of these 
f^vatloiie aitk an appropriate esef^iter program* The stmt *f>aracts*4#eics 
of Uip mqmtim* soold ikm be deteretaed f rc* a gurir of the f ri^msU ee 
of t^eir mi, IVlLs is dlwrusaM in ^«ti*s 7 ai a project reoca^^t 
f-r /dbds at 

T5i» fov^gpiri | > e^a|^9^a e®«*«iwed with lb* 1 — mil r*j — rf 
Irtaitto mar -irithoct ref**** f ,Tp t# e f e^i that st i ues in _n. easi jv # 

Is (f ) *111 itwi iqfctsqeMSt errur*. |» eeBJffrss thj * iTf*st a 
dlff*r^\ isv^-'pjh <tli ►* 



f 






in tw Vi* mueLtm (G), 

5 ■ *.,«■'[;'** > - *t r -> .] ♦ r '*' ** ' < » • 

f, U * Jr % . V*1*J» J rcjyu*. rf t*. initial Hlwi -rnnUe •* r. 

’ ( *. • 1 • Vk. Lrm. eti<w. «rr*ir j**Lr .< r* k# 

Trrr. - •' , l‘l) U t%# rm wdrff sew. 

V* -Wb* ©a,i) * 1(h) * J(i). 



Nfira». * 1 * y, - y, , •« *-ttr*ct (9) Cr*m ( 1^* > i».j tut*d* 

*», )*,{((, '-r. )]*e(k r ) 

t£at f.Cx.jyJis dlif«rentiable with respect to y Li 
tfi* interval (y, ,y^ ) w- mg< us* the «an value theorem to write » 

(*• ; e,^r*,,o. 

/,) does not vaiy greatly we may let t, (x ) ■ K and writ* 

4 '* K, -6i , + «(*•»■»)• 

Collecting trots 0~rk KV) t, r U+ £ KV)*,., - (l- i* V) 4,., * £Y„ 1) 
or, fia-h'*)*, sfa** - kV)£.,_, - + 

In the solution of the initial value problem (1), K » ^(x,,!*) » -1, 

Bence, (12+h l )6,e (24-lOh 1 )^- (la+h 1 )^.,* E(h,l). ( 15 ) 

Under the aa-.-5.tl0n that E(h,i) Is relatively insensitive with 

respect to i, we may proceed as we did in solving (9). (A# applied to 

this oase the truncation error T(h) « Ch*y *((),) , and varies 

•lowly when the h increments are mall). We will henceforth 

K(h,i) to be a o«.stanti Eh*. 

We acst solve (12 + h 1 )^ « 24-lOh 1 )«, - (12*h* )«< ,+ Eh* 

«• » 0 (no error initially), mB 
*>< • f (assumed to be non-cero) 

The solution consists of the ^*>eral solution o. 

hemigeneous part ■ 0. r7+ 0, ^ ml a particular 




e 



41*9 t* til* I* 

r ‘-^1 in u«»J/c Qr.*t'. in * ko«UH *, 

u>Vu * cumd* 1 )* - aanfc* u * 

£><* - 



>• i«r»i*Oar «elii!dft 
Uuil 



: ( r, *i J ) ( 04 - ^ 




- oi s 
ro , . ’ ' * i 



tn® w^le^e elation i« then o r % c> r[ +5 h . (U») 

Let r ■ a + it, r^» a - lb* 

Denote c = <1 + i£, - i£, where ^ h* end 4 : jV 4 -*>JL K* [ 

The solution (16) is transformed into trigyec**tric form hr Lett** 



0 * arc sin 



Jl. 



tT^TT 5 



Than, € « |rl (2*oos n0 - 2£s±n n9) + XI h* ♦ Kspandingi 



4* = f r/* h V oos n0 +7 ( /+ (a-l)~ h^) sin n0 + ~ h 
In particular, if we approximate 

a = 1 b * h ^ i.e.^r, « 1 + ih 9 r A * 1 - ih 
4* * (fl+K 1 ) ^ ^ h oos nO + ~- ( / sin n0)~V h 1 



where 0 =* arc sin 



K 



For small h 



(17) 

<«> 



f l+h 1 ) — h^oos nh +^sin nh"? + — h ♦ 

( *1 h ) 1 * 

This solution of the error propagation equation qualitatively 
plain* the behavior of the eocperip^tal res tit a displayed in Fig. 16, 

a. If h is small, the sine tea of i - Jl h” oos nh 4~fin iw*7 

C 10 K ) 

1 cLgnificant an^ the oodi»e t ere is negligible. pTr»r would be 
expected to vary nearly sinu^idally. 

b. it* ooid b- atT^.iTcijfctely equal to ter. 
ti^le -aT^, dL^Oe tb* ten-£ *L* <=* will vi*idi at Umw f«&Ms, Jivee 
h * •* - ^lHjt f-r n • if , k • t, a. .... 5 , w» 1 — T|n i «e t ta ttaa 



k 

•A 1* 
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4H4 -r**r*bmm a a ij>| All J8H Ftf> i* Awi 

to IH# C4p« S'hff iiiHlS £M la fef®4 33 «*«. *»1 ^ -«1 

cWiarty be Tor n *ttb h * .01 UJ .Q;* f ree* — 

c. “> brin Ireted Ur* In (lA) ie <yirl /*-?*> 

tfco tr»: t*rta stem with lxuvvMidac n if ill / l «rf ilailA 

lffrkl, Ttm fc-clifioatAw factor (/fS* ) in (10) it ml/ as - 

motion of J r i T Hie actual roots were foesdi to be 

. j, - ( II- fly ) 2 / Vl" 4" 4- ~ * h * 

/ s^l * C I 4 *- - * ^ 'ISh * ") 4 K *( jOa* - ^ a ^ 

!♦* * 

for all h* ^«rev«r h is snail enough to insura that (18) closely Af**-**- 
lmatas (1?) 9 arror marina should r amain nearly a distant* fig* 33 *i**» 
•tantiataa this f*r h = .01, Increasing arror marlma in Fi(e* 19 *1 

(h » *02 and ,05 f raspaotivaly) ara attributed to the failure of h to neat 
the criterion stated above* 

Vo explanation is offered for the high frequency error oeolUatioi 
so prevail t in amna of the graphs 9 a.g* Fig. 16* Investigation df this 
T#!qr*mH*iep is aruggmibed as a project for future stud^y* 

4. Pa script of the Ter? Program . 

This >n<n« it deQgaed to solve a given eeoond order differen- 
tial aquation 3 ti^# 9 spaces lively using each of 26 Pad£ o*rtr ergots 
to |pn * rate a central recursion formula* Integration of the glean 
•qua!!*!* i* the pr»U c tor-corrector method, using the SGlfte pr*rii*lrr 

forwJ.j 1 * *££r>l vt|.^ J . »*C..S*I T 

• »trt4r/ T*lu.fc nr* c— -<>'•••' V»r <Vr JU=»# .-K ,»tta — ***1. f>. lt*» -cAt. 
P?ee«i» wrr fn cirmi -at -- a- ,»*^r »• 900 ‘ i on .irar l*&v 
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^♦!*r*lr*_* at *mm Lf — H l>- f**™* -*«i s4 »* 4* 

Known , I lit plotted With %TTX>r t * *r-f!4.*wT»t *■* iilL« S< jfeiV^ «•! 

as abscises. riff. 3 miplalni ts« gp^els MKi Is Ite fflp^i* tl|j I 
is a printout of the proffr*-» to evaliate t*»« gersMUag 

tha ooatne ft»tion, y* + y * 0 {Frugfrvm HEn 2). Fig. 5 Is a f r£^w*\ 
of tha program oaad to evaluate tha Coulimf* Wat a ^imUsa 
(P ro<Y*» TMCD 6). 

5. DireeaeiOtt of Data . 

a. Program TMCD 2 

(1) Pigs. 7 through 25 are graphs of progressiva « rr*\ - limit- 
ing when tha initial value problsm y" + y « 0 t y(0) * l t y f (0) > 0 U 
solved using an integration step length of .02 and corrector tmvmC u 
derived from tha various Pad£ coovargsnts. Pig. 6 relates tha r a u^ei ^i 
formula number shown on eaoh graph to a description of the oorrs^ee^ri 
Pads fraction. Formulas l6 v 17, 20, 21, 23, 24 and 26 diverged regldly 
and no associated graphs were plotted. 

(a) Formulas 10, 11, 12, 15, 22 and 25 produced remarkable 
accuracy over the entire range of 900 iterations. 

(b) Formula 6 (Pig. 12), although less aomsrate ihm 

group Just dieiussed, has an ijf^rtant property: error »***»* are 

diminished with increasing iteration. 

(o) The remaining graphs indioate that error are 

not being da»-*d yet are increasing linearly with regfmat to n, tbe 
n*mtmr of itmreti Jli, 

(2) P**e sal cf data was obtained by repeating the nffm 

af t*r d*orvL*«n# th* ievweti^ s*m lemffOi to .91. 
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Th* tollo-tni 

IS# *MM P"0*«3 !t ft* majt ^ «r*C* #tW» • *t^ 

ler^th t»f .EC »•*•**** It}, li, 11, i;, S ••>-* -*'■'• f--«« w t • 

l ... •teiU'.; « wit r •->.» mAX*<- iwttt w> vat u* 

«a« ned). f*mli« 1)| t? «*J 21 sbcv t «Mll *r^*T •*JUTl**tirw v*i| • 
til U)« r V. -ml* <T4V e I»h crater srrcr v^UTl'^Uaa. 

fUi, ?P tnrr J> j) *r» «*<-!'< tfc*ee rw*2t* 

(b) r • r ni* 6 pecwlrles L#«e *oo--r »t* result t tt* nit *M* 
#re*t*- <>t»p lerc** bst it* error dicing oh*rect*ri#Uas m*.ir (Fig. ^7) 

(c) Formil* 5 is «>rw *sefcle. Its sooorsoy is j'Wwtll* 

to thet In the preceding p*. (Fi«. 26) 

(d) All other f orrwiles proved to be erre soeurste Ub4 
u*ed with the gre-ter step x— - th. All of thst, however, saplify error*. 

(e) The • rror prsWbs a*ooeiated with fonaUi 9» 10, 11, 

13, 14, 18, 22 aad ?5 T T-re i<J*otioel. (Fig. 33) 

(3) Increeeirr the step l^a^th to .05 hed th* following effeotst 

(s' Th’ forvules whisfc proved best with s step lwsftb of 
.02 ere lee* eoourete with • step length of .05. Figs, y* thmgb 
41 sr* th* 1 error jr-ephs as*<r*:l«ted vita Foi«iAae 6, 10, 11, 12, 13, 15* 

22 end 25 <wtVr* the .05 «Mg> l-«*glS*. The*e gra**s «*ow the only fo*«iles 
rwlairvln* *ary error d—nlng o(*aracteri«tics to be mtfers 11, 12 gpl 25. 
b. P*v<r*A. 2JCD 6 

((j ftw* Kste (|llfQ4B •»» snivel for taw c<m*j1£*.*1**h 

if d -at v l .VyOl.-idi t*.V 1 2, y.J)* *Jt* ••rtc-x# le*1 e- 

* ef> i teethe. ft •*« 0*1 feealM* 1* fWtSWU tfNT n*wb* fee 
*iM - n t#. r ’e. 42 (let* #**?*-> r«2e Sf ivwUUor. -J zs-puimi *&1. 
itt’S* C *:• e>* Ot ess ifltf v— -i «.r *f lUw»U«* istlctltJ. Im* 






I 















i? 



Isiuti 'to* -JK ai wr» « M » w» 

l»- (l 1 10 1 ). UA.«t’iUi» ftf fir. <M itoi !aUMU|i 

(«) 1*«14 o 41 r«4ults awe abt*Ui«d Lj. *,U **.►-- 

f«irMi.ju 9 , 13 , 14, 13 «ru 18 w«r« u«J. tW * £n>.* tlr*- ".£» 

Mit fccc'ir tt* solutions. :'sd*c-in< rtsp Lj*>tot< r 'iur»ij \Vrc^>. 

out the range. 



(b) Formula 8 solutions wrte slightly less **»o*Lr«,te Ltan 

tho«e of the group above in solving the equation y" + y (1-i-i ) * 0, 

« «» 

Solutions of the squition y" + y( 1-- -1) = 0, using s step lengta 

* * 

of .02 were more accurate than with any other formula. 



6. Conclusion*. 



In his Numerical Methods [Yj f Kopal Illustrated the shortcomings of 
going-ahead recursion formulas: those derived from the relationship 

f K 7 1 

y^ = ^ . When using going -across formulas stemming from the 

Integration formula y t » )~y~~ | 8111 of f, ®st be ob- 

tained before y t can be calculated. Therefore, the method suggested by 
kopal is essentially the predictor-corrector method ooupled with a pro- 
cedure for deriving an unlimited supply of corrector formulas. Some of 
these we have seen produce quite accurate solutions to several differen- 
tial equations of second order. A glance at Figs. 43 - 45 f the error 
of y" + y a 0 after eolation by the Milne method, is evidence of the 



lority *t a ntefemr of the correctors derived from Psde conversant#. 

It is *«^eeted at this point that the les* accurate gtlnf -ahead 



tyre fr^ulae ni 4 x»t be generated in the tew way and ^sed as predictors. 
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fat (Mt 



• 'Wt* 



m «ll ftrpridv ■4Tl*4*rt k-Mnart <** * *4*4 1»* 

«j<Q[ at ►rail tC.r 4»ni«Hai, 

vai# «C t/» *ott*o*cs !«*qui <1w i Ui Til* jMp* a* *u*u 

V*t *lj *,-4ttcnr k, ' :rufw. of fu KkUUaMl fc- m)i> o—d< I- 

fn» i ptSjnaA*] Ul| t *f 1 1 it r 

t— a- — • i«w nut aty-*i . In 4*aU»C wilt **■**) atorf 

oMipaWrl stttf Kl|j**r 4 «yr *i It will b* srT.SMwV 

to Ioo*t* val **» rf p*r*-*.*ra rd«4 Bill ratoA* *11 root* MUm th* 

0-* «s**ll«t raSktlt* obtained with Ffra&s 15 (FI*. W) V^±5 4*fc* 
that * fn^rth m&m pnijrea«L&l innojetia to £ ^ i* mt+ynm'i 

to u>au*% to*t truncation **ror i* t*»t * «if.fi**nt ftoM**. tv wrwll 
*-»”•■, 'Jtarafor*, that la* tmm of polyno*l*l« of hisaar d-er** ***,14 **— 

♦ It* fertdltltcwl o»«ll»*tl«v*. 

In *ov*T 7 , it Is csncludad that th* prooadur* saffastad ny a**akl 
sfwnid b* «*») at * aatbod sf terl*4j)( pradlotor and on wi m fiUi 
2b* f*«lc* *ahi***«l wiio so** tS th* corraobors Oarirad 1* this *aW 
lads oat* tfiat ttd* procedure i* (rodstjif Asd Merits f>irto*r tUnts 

?. fc- /X JLafc C--Jt-'r'±i* 

• . Dw**tt*ata til* roots af tl* *X*raeb*ri*loo assist*! 

■sist tbs ass fi-r*4 dat reeve- also frnail n rH arms ml ia ibis la* 

"•ruix «od t*be*abs Tal'*aa a£ I* far Shi** all root* Ua within th* «Ut 
hlti, tbit *t* . **11 >01*, .st* »ctM raadll^ s*i**t to* boat 
l*> l'M« fto*«l» far sol*!/* * l»*> afdlUai ostr a 







r«<r* < f o< ** * -tM •:<«, 

• 

... ZrnfU, . ,‘w« t*bJ. nbwa ■W lft « . * :• ,' - '**• * i 

!r* m n,\*zh m *v-*in * ttf |%X^ iNii rmiiA 

ptj miIm- AiMLh r u*T* to tmn tmt tHi I^u . — i*» • 

C*4c._Lete •rniT* «»*■» wit^» >tq, 4J, 0* fVI * 

■ ■a# »rlt* ut r “Victor •’ *ri?e‘ [m rMixi) i * <* 

ctx.7*tx9m ca. lalt^ff] duprlbe* ti»o edVAAfeag** of <s&fq pmAl staff* 
fif'd corrector! of the wt order. 

( #. « , ^ e* i * 

c. Extuid the polynomial repre^ntetion of 5 — 1 f^M*! 

. » K * 

term* esid derive higher order recur^on forrsiii#. hp e*tAf 4 t« 
the roots of these formulas as s^Qiested 3n subpArefrtp* a a!«rre* 
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# » f%V* 41 ll'#> »*T ifl* Stale ^ u»lt «|Hl&l^ 

h. &*n«» *. P«d» IMA* wU- a »*I*H/W»»* tmt "• '* 



On* <#*1 ca i - » mm mt *f r>W-A#*J r^ Ane* fwW, W*. 

ttlLtiU «•-.. of If. f»rt far tfe© fr-Hit r» 1 .*-«*!* um lr* WS’% 

T «3T 6. C41frAU<* *rrofi «M a:-x»kr« wi»A ft<. fc£, 3>. 'All *«l<m . 
••*l'j»r-* **.. -tmrt is of (V^fclictor forvuM f i ••» I 

o.rv«r,-as. ♦t-r, J" } d«*crl> • i t*>« *i <«nWr»« '~ f 'Mm p»»41c 4 rfi 



*nj correctors of th® •'■r® ord#r. 

o. Ext *ed the polyn:-dal reprer« tat loti of { 

•wen t»c&m and derive higher order recursion fomtiLes* 





the roots of these formulas as pif tested in subpar^regft a 
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•v . ) * L* --f VT^sr^ ^1* 
c ttriiMo 



D«*r--r^^t variable 



rF 

rop 

TT 

2?l** 

D(l) 

Z(I) 

AN(2) 

AD(4) 

^X, AI> 
ATP, ATDP 

JP 

JQ 

H 

N 

M 



y f 

y" 

y 

y - y 

▼y 

▼V 

Coefficient of 2 nd nweretor term Ln Fade c+ jgv w y asit 
Coefficient of 4 th denominator ter® ln Pade *ot?*«skK 

Initial values 

Number of first convergent to be evaluated (Hjp*fl£) 
Number of last convergent to be evaluated 
h (Integration step length) 

I? iber of starting values to ta computed ty W>£e~iiitta 
Total mjeber of iterations desired 



EN(I) Denotes number of terms in n*»rator of foMJJa mo. (1). 

Used for information in printout 

LD(I) Smet as IN(I), except refers to dan^inatar tar « 

!KL(K t <* y J) Array which stores Pads coefficients, »■ 1 donates 

naaerator coefficient; K* 2 , daft'jctsator coef?ici**t, 

I is index no, of convery^t. J is the n*^mr oi a 
K*ecific t«ri # t,g, # HL(2 # 7,4) dp*^ota« T -ha ^n#Uec 
of the fourth ter* in the d*r*«Uai*r a* £^# 7 ^,! >*.7 



QO Coefficient if *, in reC'«Aes fUfflgl ■ 

tf(j) Cofcficieot of #, ; in rwrdcr, focwel* 
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